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ABSTRACT 


A  simplified  version  of  the  shock  expansion  theory  is 
employed  to  estimate  the  unsteady  aerodynamic  load¬ 
ing  on  a  planar  surface.  Expressions  for  oscillatory 
aerodynamic  coefficients  are  developed  for  the  case 
of  a  three-degree-of-freedom  airfoil  section. 
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Free- stream  density 
Shock  wave  angle 
Circular  frequency 


Conditions 

Conditions 

Conditions 

Conditions 

Conditions 


immediately  downstream  of  shock  wave  at  x  = 
immediately  upstream  of  shock  wave  at  x  =  x^ 
on  the  lower  surface  of  airfoil 
at  the  leading  edge 
on  the  upper  surface  of  airfoil 


I.  INTRODUCTION 


The  two-dimensional  shock  expansion  theory  has  been  employed  to 
estimate  the  steady  airloads  on  slender  airfoils  traveling  at  h^^rsonic 
speeds  (see,  for  example.  Ref.  1).  An  unsteady  version  of  this  theory  has 
also  been  proposed  for  the  determination  of  transient  and  oscillatory  airloads, 
on  slender  airfoils,  in  the  hypersonic  regime. 

The  unsteady  theory  may  be  developed  from  the  hypersonic  small  dis¬ 
turbance  equations  for  unsteady  flows  (see  Appendix  A)  or,  if  physical 
arguments  are  employed,  directly  from  the  results  of  the  steady  flow  theory. 
The  resulting  formulae  are  rather  complicated,  but  they  may  be  simplified 
when  the  time -dependent  disturbances  are  small  compared  with  the  steauly 
disturbances.^ 

The  procedure  developed  in  Ref.  Z  is  briefly  outlined  in  the  present 
work,  and  expressions  are  developed  for  the  pressure  distribution  in  various 
flow  regimes  of  technical  interest.  These  results  then  are  employed  to 
determine  aerodynamic  coefficients  for  a  two-dimensional  airfoil  oscillating 
in  the  three  degrees  of  freedom  of  plunging,  pitching,  and  control  surface 
rotation.  Some  example  calculations  are  carried  out  for  a  flat-plate  airfoil. 
The  assumptions  upon  which  the  theory  is  based  are  noted  in  Appendix  A. 

The  sign  convention  employed  in  this  paper  is  as  follows.  The  flutter 
sign  convention  is  used  in  the  unsteady  case:  unsteady  forces  and  deflections 
are  positive  down;  rotations  are  positive  with  the  leading  edge  up.  The  aero¬ 
dynamic  sign  convention  is  used  in  the  steady  case:  steady  forces  and 
deflections  are  positive  up;  rotations  are  positive  with  leading  edge  up. 
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n.  AERODYNAMIC  THEORY 


A  simplified  version  of  the  shock  expansion  theory,  developed  in 
Ref.  2,  is  employed  to  determine  pressure  distributions  inflow  regimes  of 
technical  interest.  Shock  expansion  flows,  expansion  flows,  and  flows 
downstream  of  secondary  shock  waves  are  treated. 

Shock  Expansion  Flow 

Cons ide  r  the  hype  r sonic  flow  ove  r  the  uppe r  surface  of  the  configuration 
shown  in  Fig.  1.  The  surface  slope  at  the  leading  edge  is  positive; therefore. 


Fig.  1.  Coordinate  system  for  slender  surface 
in  a  hypersonic  flow. 


the  flow  is  compressed  through  an  oblique  shock  wave.  Conditions  are 
assumed  to  be  such  that  the  shock  wave  is  attached.  Assuming  no  appreciable 
compressive  turning  of  the  flow  downstream  of  this  point,  we  find  the  approxi¬ 
mate  formula  for  the  aerodynamic  pressure  on  the  upper  surface  to  be 


p(x,t)  =  Pjj(t*)  1  +  M^(t*)[  6(x,  t)  -  6^(t#)lj 


2y/(’y-l) 


(1) 
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where  6(x,  t)  ia  the  flow  deflection  nngle  at  point  x  on  the  surface  at  the 
time  t.  The  qiiantities  and  M^Ct*)  denote,  respectively,  aerodynamic 

pressure  and  Mach  number  immediately  behind  the  leading  edge  shock  wave, 
evaluated  at  the  time  t*  s  t  -  x/J3^.  The  flow  deflection  angle  at  the  leading 
edge  is  denoted  by  6^.  The  pressure  Mach  number  M^(t^)  depend 

upon  the  flow  deflection  angle  at  the  leading  edge,  at  time  t*.  and  the  free 
stream  conditions.  They  may  be  evaluated  from  the  following  formulae: 


P«  =  +  VP^ki  6 

n  oo  00  oo  n 


^oo®n  +  [  ^  *^oo®n)  ] 


211/21 


K-j 


(y  +  1)^M^0^ 


2yM^02  -  (^  -  1)]  [(Y  -  DM^^^  +'^ 


oo' 


2ll/2 


.  (2a) 


(2b) 


(2c) 


where  <f>  denotes  the  angle  that  the  leading  edge  shock  wave  makes  with  the 
free-stream  direction.  These  expressions  are  the  approximate  form  of  the 
oblique  shock  wave  relations,  and  they  are  developed  under  the  assumption 
that 


2  2 

sin  6  w  6  ,  sin  p »  p  ,  6  «  1  ,  M.  »  2 

n  n  n  oo 

Specializing  to  the  case  of  small  unsteady  motions,  we  can  simplify 
Eq.  (1)  further.  As  in  Ref.  2,  the  flow  deflection  angle  6(x,t)  is  written 

6(  X,  t)  =  9(x)  +  «(X,  t)  ,  with  I  €  I  «  1 0  I  , 
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where  c(x,  t)  denotes  the  time -dependent  contribution.  At  the  leading  edgOt 
the  flow  deflection  angle  is 


6  s  e  -f  c 
n  n  n 


The  aerodynamic  pressure  and  Mach  number  at  the  leading  edge  are  esqpanded 
about  their  respective  steady  values  p^  and  Retaining  first-order  terms 

in  these  expansions  gives 


pj^(t)  =  P^[  1  * 


M  (t)  =  51  [  1  +  m€  (t)]  , 

n  n  & 


where  t)  and  m  are  the  dimensionless  rates  of  change  of  pressure  and  Mach 
number  with  turning  angle. 

The  expression  for  m  are  obtained  from  the  oblique  shock  wave 

relations: 


2y(y  +  1)M^(M^^' 


[2y(M^?)^  -  (y  -  1)]  +  l| 

n  <y*  ■>(“»») 

^  ((V  -  +  1 


Equations  (3)  are  now  substituted  into  Eq.  (1)  snd  terms  of  t  are  neglected 
so  that 


p{x,t)  =  pjl  +  +  «(x,t)  +  €^(t*)[mA(x)  -  l]|J 


Zy/Cy-i) 
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where  ^x)  =  6(x)  -  6^.  Equation  {5)  may  be  expanded  if  we  aasume  further 
that 


A(x)  + - - -  »  €(x,t)  +  «  (t*)[mA(x)  -  1] 

(Y  -  ”  ' 

Neglecting  terms  of  and  higher  in  this  expansion,  we  obtain  the  following 
relationship: 


p(x.t)  =  Pj^[Fj{x)]2v/<Y-1) 


1  + 


Titbit*)  +  |y-^j|*(x,t)  +  €^(t*)[md(x)  -  l]| 


(6) 


where 


Fj(x)  =  1  + 

when  Eq.  (6)  is  split  into  time -dependent  and  time -independent  parts 
p'(x, t)  and  p(x),  the  respective  contributions  are  seen  to  be 

p(x)  ,  {7a) 

p'(x,t)  =  iR  j(x)  €(x,t)  +  31 


where 

3{j(x)  = 


(7c) 

(7d) 
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The  quantities  and  are  obtained  from  Eqs.  (2)  by  replacing  6^^ 

with  0  . 
n 

Expansion  Flow 

In  the  case  that  the  surface  slope  at  the  leading  edge  is  negative!  an 
expansion  immediately  takes  place.  Assuming  no  appreciable  compressive 
turning  of  the  flow  downstream  of  the  leading  edge,  we  find  that  the  approxi¬ 
mate  formula  for  the  aerodynamic  pressure  on  the  upper  surface  is 

P(x.t)  ^  Poo[l  + 

The  flow  deflection  angle  is  written  as 

6(x,t)  =  9(x)  +  €(x,t)  .  with  |€|«|e| 

and  it  is  assumed  that 


1  +X^M  e(x)  » 
2  00 


Equation  (8)  is  expanded  and  terms  of  order  €  and  higher  neglected  so  that 


P(x,t) 


1  + 


vM^€(x,  t) 


T^xT 


(9) 


where 


F2(x)  =  1  +  X^M^e(x) 
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Equation  (9)  may  ba  split  into  time  •dependent  and  time -independent  parts 
p'(x,t)  and  respectivelyt  where 


p(x)  =p^[r2(x)]*^'<V“«  ,  (lOi,) 

p'(x,t)  =*3(x)  €(x.t)  .  (10b) 


and  where 


»3(x)  =  (11) 

Downstream  Flows 

Modification  of  the  previous  axuilysis  is  necessary  if  the  unsteady  aero¬ 
dynamic  pressures  downstream  of  a  second  shock  wave  are  required.  Such 
a  situation  arises,  for  example,  in  the  determination  of  the  aerodynamic 
pressures  on  the  compression  side  of  a  deflected  control  surface.  The 
difference  in  the  analysis  arises  from  the  fact  that  the  flow  conditions 
upstream  of  the  shock  wave  are  now  unsteady. 

The  expression  for  the  surface  aerodynamic  pressure  in  such  a  flow  is 
derived  under  the  assumption  that  the  surface  downstream  of  the  shock  wave 
has  negligible  curvature.  Two  types  of  upstream  flow  are  considered, 
namely,  a  shock  expansion  flow  and  an  expansion  flow.  The  analysis  pro¬ 
ceeds  under  the  assumption  of  small  unsteady  motions,  and,  since  the 
development  is  similar  to  that  of  the  preceding  sections,  only  the  major 
results  are  presented. 

Consider  that  the  oblique  shock  wave  originates  from  the  point  3^  on 

the  upper  surface.  The  aerodynamic  pressure  p  and  Mach  number  M 

c  c 

immediately  behind  this  shock  wave  are  given  by  Eqs.  (2)  with  p^,  M^,  p, 
and  6^  replaced  by  p^^,  M^,  and  5^^,  respectively;  where  P|^  and  M^^ 
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denote  the  aerodynamic  preeeure  and  Mach  number  immediately  V9»ttm§m 
of  the  shock  wave  at  Xj^,  and  where  and  denote  the  shock  wave  angle 
and  flow  deflection  angle,  respectively,  both  angles  being  meaeured  from 
the  flow  direction  immediately  upstream  o£  Om  shock  wave.  "Rie  pressure 
distribution  farther  downstream  of  fee  shock  wave  is  determinod  from  slmck 
expansion  theory.  These  results  may  be  simplified  for  the  case  of  small 
unsteady  motions  by  expanding  the  various  quantities  about  their  respective 
steady  flow  values.  When  this  procedure  is  carried  through,  the  following 
expression  is  obtained  for  the  aerodynamic  pressure  on  the  surface  down¬ 
stream  of  the  shock  wave: 

p{x,t)  =  P(X)  +  p'(x.t)  ,  (X  >  Xj^)  , 


where 


^x)  =  .  (12*) 

p'{x,t)  =  VPh^CjF3(^)  +  CoPf^it^^)  +  yPi,W^Co[«(x.t)  -  .^(t^)]  .  (12b) 


and  where 


Cq  =  1  +  +  *'4) 
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The  barred  quantitiee  denote  steady  flow  values.  The  unsteady  aerodynamic 
pressure  immediately  upstream  of  the  shock  wave  at  is  denoted  by  pj^. 
and  the  unsteady  contribution  to  the  flow  deflection  angle  immediately  down¬ 
stream  of  this  shock  wave  is  denoted  by  t  .  In  the  above  formulae  both  of 

c 

these  quantities  are  evaluated  at  the  time  tj^  =  t  -  (x  - 

The  quantity  v  is  determined  from  the  first-order  expansion  of  the 
Mach  number  about  the  steady  flow  value 

Mj^(t)  =  1  +  v{t)] 


The  various  terms  appearing  in  Eq.  (12)  are  dependent  upon  the  flow  condi¬ 
tions  upstream  of  the  shock  wave  at  x^.  If  these  flows  are  specified,  then 
the  equations  can  be  reduced  to  a  more  explicit  form.  With  a  view  to 
technical  applications,  two  particular  upstream  flows  are  considered,  namely, 
a  shock  expansion  flow  and  an  expansion  flow. 

Shock  Expansion  Flow  Upstream 

In  this  particular  case,  the  analysis  of  the  first  portion  of  Section  II 
can  be  employed  to  determine  p^  and  v(t).  Substituting  these  results  into 
Eq.  (12),  we  find  the  unsteady  aerodynamic  pressure  to  be 

p'(x,t)  =  3l|^«(x,t)  +  *7{*h)‘n^**^  '  ^  *h^ 

(13) 

where 

«5  =  -3I4  ,  (14b) 
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(14c) 


I  nr  3^ 

1  Y. -J. 


S^(x)  =  Cq»2(x)  + 


m 


^_jITJl-mA(*)l 

+  ■*-2 - rjsi — 


.  (14d) 


The  coefficients  91  ^  and  9(  ^  are  evaluated  immediately  upstream  of  the 
shock  wave  at  Xj^. 

The  coefficients  fR  (x)  and  91  -(x)  are  defined  by  Eqs.  (7c)  and  (7d), 
respectively.  The  steady  ilow  aerodynamic  pressure  ^  and  Mach  number 
immediately  upstream  of  the  shock  wave  are 

The  steady  flow  Mach  number  immediately  downstream  of  the  shock  wave 
is  determined  using  the  appropriate  form  of  Eq.  (2b).  By  means  of  the  above 
equations,  the  unsteady  aerodynamic  pressure  downstream  of  the  second 
shock  wave  can  be  related  to  the  undisturbed  flow  conditions  and  the  surface 
geometry. 

Expansion  Flow  Upstream 

In  this  particular  case,  the  unsteady  aerodynamic  pressure  is  found 

to  be 


p'(x,t)  =  »:g€(x,t)  +  9ig€^(t|^)  + » ^  *h^ 

(16) 


-11- 


where 


*'8  = 

*•9  =  VPij^Cj  -  *3  ,  (17b) 


« jo(x)  =  »l3(x)Co  - 


(17c) 


The  coefficient  91  is  evaluated  immediately  upstream  of  the  shock  wave  at 
and  the  coefficient  iR  3  is  defined  by  Eq.  (11).  The  aerodynamic  pressure 
and  Mach  number  are  given  as 


Ph  =  Po 


-\12v/(y-1) 


(18») 


^=^oohK)r' 


(18b) 
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m.  PRESSURE  DISTRIBUTION  ON  AN  AIRFOIL  STRIP 

TTie  results  of  Section  II  are  now  employed  to  estimate  the  unsteady 
aerodynamic  loading  on  a  two-dimensional  airfoil  strip.  (Such  an  airfoil  strip 
is  shown  in  Fig.  3a.)  The  leading  edge  of  the  airfoil  is  located  at  x  >  0,  and 
the  control  surface  hinge  line,  at  x  >  Xj^.  The  various  loading  cases  have  to 
be  distinguished  depending  upon  the  conditions  at  the  airfoil,  nose  and  at  the 
control  surface. 

Case  A 


Case  B 


Case  C 


Case  D 


The  additional  subscripts  u  and  i  refer  to  the  upper  and  lower  surfaces 
of  the  airfoil,  respectively.  Shock  waves  are  assumed  to  originate  only  from 
the  leading  edge  and  from  the  control  surface  hinge  line.  Ebcpressions  for  the 
net  unsteady  aerodynamic  pressure  acting  upon  the  airfoil  are  listed  for  these 
cases.  These  expressions  are  modified  for  the  case  of  harmonic  oscillations 
and  are  simplified  in  the  case  that  the  reduced  frequency  of  oscillation  is 
small. 

If  w(x,t)  denotes  the  unsteady  lateral  motion  of  the  airfoil  strip,  than 
«y(x,t)  =  -s^(x,t)  =  -^-^  + ipJ- ,  {OSxSZb) 


The  four  particular  cases  to  be  treated  here  are: 


Q  >0 
nu 


Control  surface  deflected  down 


®nl>0 


0  >0 
nu 


Control  surface  deflected  up 


«  ,  >  0 
ni 


0  <0 
nu 


Control  surface  deflected  down 


®ni>0 


0  <0 
nu 


Control  surface  deflected  up 


®ni>0 
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The  net  unsteady  aerodynamic  pressure  p^Cx^t)  is  defined  as 

p»(x, t)  =  p|^(x, t)  -  p'n(x,t) 

Small  Unsteady  Moti<ms 

The  net  unsteady  pressure  between  the  leading  edge  and  the  control 
surface  hinge  line  is 

^ —  =  -G^t  I  (x,  t)  -  G^€  |(0,  t*)  ,  (0  <  X  <  Xj^)  I  (19) 

"oo 

where  m  =  li  n  =  2  for  Cases  A  and  B;  m  =  11,  n  =  12  for  Cases  C  and  D. 

The  net  unsteady  pressure  between  the  control  surface  hinge  line  and 
the  trailing  edge  of  the  airfoil  is 

.  -c  (X,  t)  -  O  (0, 1.)  -  tj-)  -  a,.,  (xj,  .j)  . 

*^00  ^ 

(Xj^  <  X  <  2b)  (20) 

where 

p  =  3,  q  *  4,  r  =  5,  3=6,  for  Case  A, 
p  =  7,  q  =  8,  r  =  9,  s  =  10,  for  Case  B, 
p  =  13,  q  =  14,  r  =  5,  s  =  6,  for  Case  C, 
p  =  15,  q  =  12,  r  =  16,  s  =  17,  for  Case  D. 

Explicit  expressions  for  the  coefficients  G^  are  given  in  Apj>endix  B.  The 
argument  x^  indicates  that  the  limiting  function  value  as  x  -*Xj^  from  above 
is  to  be  employed,  and  the  argument  xj^  indicates  that  the  limiting  function 
value  as  x  ••x.  from  below  is  to  be  employed. 
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Small  Harmonic  Motions 


When  the  airfoil  strip  performs  harmonic  oscillatims  of  frequency  ut 
the  unsteady  displacement  may  be  written 

w(x,  t)  B  W(x)exp(iwt) 

Therefore 

€  Jx.t)  =  -•jj{x,t)  B  -T(x)exp(iwt)  , 


where 


•n  \  ^ 


k  = 

k  -  -p- 


and  where  b  denotes  the  semichord  of  the  airfoil  strip. 

The  expressions  for  p*(x,  t)  corresponding  to  the  case  of  harmonic 
oscillations  may  be  obtained  directly  from  Eqs.  (19)  and  (20)  by  direct 
substitution  for  <  | . 

Low-Frequency  Oscillations 

If  the  reduced  frequency  of  oscillation  k  is  sufficiently  small,  the  ' 
exponential  terms  of  argument  -ikx/b  and  -ik(x  -  X|j)/b  which  appear  in  the 
expressions  for  p*(x,t)  may  be  approximated  by  power  series  expansions. 

If  terms  of  order  k  and  higher  are  neglected,  then  the  expressions  for  the  net 
unsteady  aerodynamic  pressure  are  found  to  be 


t)  *  £«>'  *  I  -  Tr(lj 


(0  <  X  <  Xjj^) 


^  •  -o/(*)  -  o,|2{21(ik  ♦  l)  t  ^  [i 

-°.|rK)  (*^  *  -^)  *  ^K)  [‘ 


ikx  (x^ 

T"  ■  tie: 


ik— B - J 


X  - 


ik 


*h  k- 


<^‘]| 

f^)1 


(Xj^  <  X  <  2b) 


(22) 


where  p*(x,  t)  =  ^exp(i«t).  The  various  loading  cases  are  distinguished  by 
the  subscripts  m,  n,  p,  q,  r,  s.  These  subscripts  are  assigned  as  before. 


-16- 


IV.  OSCILLATORY  AERODYNAMIC  COEFFICIENTS 


When  the  airfoil  strip  has  specified  degrees  of  freedom,  it  is  possible 
to  define  aerodynamic  coefficients.  In  the  case  to  be  considered,  the  airfoil 
is  assumed  to  have  a  rigid  chord  and  a  rigid  control  surface.  The  airfoil 
motion  may,  therefore,  be  described  by: 

1.  The  vertical  displacement  of  the  leading  edge, 

2.  The  rotation  of  the  airfoil  strip  about  the  leading  edge, 

3.  The  rotation  of  the  control  surface  about  the  hinge  line 
(relative  to  the  airfoil  chord). 

Considering  harmonic  oscillations,  we  obtain 

w(x, t)  =  w(x)exp(iut)  =  [h  +  ox  +  P(x  -  -  Xj^)] exp(iwt)  ,  (23) 

where  1(x)  denotes  the  unit  step  function. 

The  oscillatory  aerodynamic  coefficients  Lj^,  L^,  L^,  Mj^,  M^, 

T,  ,  T  ,  and  Ta  are  defined  as  follows: 
ho  p 

L  =  4  p^w^b^  I  +  *  (24a) 

M  =  4  p^w^b^  I  +  *  ^pP) •  (24b) 

where  L,  M,  T  denote  lift,  moment  about  the  leading  edge ,  and  moment 
about  the  hinge  line,  respectively,  generated  on  an  airfoil  strip  of  unit  width. 
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These  quantities  are  related  to  the  aerodynamic  pressure  t)  by  the 
following  equations: 


fZb 

L  =  I  p*  dx 
•'0 


(25a) 


rZb 

M  =  I  xp’!'  dx 
•'0 


(25b) 


rZb 

T  =  I  (x  -  dx 


(25c) 


The  oscillatory  aerodynamic  coefficients  are  determined  when  the  appropriate 
expression  for  p*  is  substituted  into  Eq-  (25),  the  resulting  equations  inte¬ 
grated  and  then  identified  with  Eq.  (24). 

For  the  remainder  of  the  analysis,  it  will  be  assumed  that  the  reduced 
frequency  of  oscillation  of  the  airfoil  strip  is  small,  so  that 


k 


2 


«  1 


The  approximate  expressions  (21)  and  (22),  therefore,  may  be  employed  to 
determine  the  aerodynamic  coefficients.  Substituting  Eq.  (23)  for  the  down- 
wash  in  these  expressions,  we  arrive  at 

p*(x,t)  =  -p^l^  Hj(m,n)  +  aH2(m,n)jexp{iwt)  ,  (0  <  x  <  Xj^)  ,  (26) 

p#(x,t)=  -p^^^  H2(p,q,  r,  s)  +  QH^(p,q,  r,  s)  +  PHg(p,  r)Jexp(iwt) 

(Xj^  <  X  <  2b)  .  (27) 


The  coefficients  are  defined  in  Appendix  B. 
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These  expressions  are  now  substituted  into  Eqs.  (25)  and  the  resulting 
equations  are  integrated.  Comparison  of  these  results  with  (24)  shows  that 


the  oscillatory  aerodynamic  coefficients  are: 

1^  =  -|i[K^(m,  n) -I- K2(p.  q>  r.  s)]  ,  (28») 

=  -p[Kj(m,n)  +  K^(p,  q,  r,  s)]  ,  (28b) 

Lp=-pK5(p,r)  ,  (28c) 

=  -p[K^(m,n)  +  Ky(p,  q,  r,  s)]  ,  (28d) 

=  -|i[Kg(m,  n)  +  Kg(p,q,  r,  s)]  ,  (28e) 

Mp  =  -pKjQ(p,r)  ,  (28f) 

Tj^  =  -p|^Ky(p,  q,  r,  s)  - K2(p,q,  r,  8)j  ,  (28g) 

=  -p^K^(p,q,  r,  s)  - -^K^(p.  q,  r,  s)j  ,  (28h) 

Tp  =  -k[Kj(,(p.  r)  -  ^  K5(p.  r)]  .  (28i) 


where 


*^oo  1 


=  yM 


2 

p 

00^00 


I 
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Table  I.  Oecillatory  aerodynamic  coefficiente  s  S) 
flat -plate  airfoil  at  eero  angle  of  attack. 


1/k 

Real  [  1^] 

Imag.  [  Ljjl 

Real  [LJ 

X]»ag.Ci.J 

Ref.  No. 

0.41667 

0.00147 

-0. 08369 

-0.  03463 

-0.08321 

3 

0 

-0. 08334 

-0.  03472 

-0. 08333 

0.  99206 

-0.0035 

-0.1965 

-0.  19843 

-0. 19557 

3 

0 

-0. 1984 

-0.  19684 

-0. 19841 

2.08333 

-0.00715 

-0.4214 

-0.  88295 

-0.40836 

3 

0 

-0.4167 

-0.  86805 

-0.41667 

4. 96032 

-0.00826 

-1.0108 

-5.01927 

-0.  97069 

3 

0 

-0.9921 

-4.  92095 

-0.  99206 

14. 8810 

-0. 00848 

-3.0370 

-45.  1986 

-2.9112 

3 

0 

-2. 9762 

-44.  2888 

-2. 9762 

1/k 

Real  [Mj^l 

Imag.  [Mj^] 

Real  [M„] 
o 

Imag.  [Mj^l 

Ref.  No. 

0.41667 

0.0027 

-0.08477 

-0.  03448 

-0. 11119 

3 

0 

-0.08333 

-0.  03472 

-0.  mil 

0.  99206 

-0.00347 

-0. 19389 

-0.  19732 

-0.26113 

3 

0 

-0.  19841 

-0. 19684 

-0.26455 

2.08333 

-0.00917 

-0.41947 

-0.  88148 

-0. 54470 

3 

0 

-0.41667 

-0.  86805 

-0.55555 

4.  96032 

-0.01093 

-1.0099 

-5.  0177 

-1.29435 

3 

0 

-0. 99206 

-4.  9210 

-1. 32275 

14.  8810 

-0.01130 

-3.03667 

-45.  197 

-3.88152 

3 

0 

-2. 9762 

-44.  289 

-3.  96826 

The  coefficients  K  are  defined  in  Appendix  B.  They  may  be  evaluated  for 
given  airfoil  profiles.  The  parameters  m,n,p,q,  r,s  identify  the  varioue 
airfoil  and  control  surface  orientations. 

A  munerical  example  will  be  carried  out  to  illustrate  tte  theory.  A 
flat-plate  airfoil  will  be  treated.  The  reduced  frequency  of  oscillation  is 
assumed  to  be  small,  and  the  low-frequency  approximation  is  employed. 

Consider  a  two-dimensional  flat-plate  airfoil  oscillating  in  two  degrees 
of  freedom,  namely,  pitching  and  plunging.  If  the  airfoil  is  at  zero  initial 
angle  of  attack,  then  the  nonlinear  aerodynamic  effects  are  absent  and  the 
aerodynamic  coefficients  derived  from  this  theory  should  agree  closely  with 

3 

those  obtained  from  the  linearized  supersonic  theory.  The  results  of  both 

theories,  for  M  v  5  and  for  various  values  of  k,  are  shown  in  Table  I. 

00 

The  agreement  between  the  two  sets  of  aerodynamic  coefficients  is  very  good. 
Ihe  small  differences  that  do  appear  arise  from  the  quasisteady  nature  of 
the  present  theory  together  with  the  low-frequency  and  high  Mach  number 
approximations  employed. 

When  the  airfoil  is  at  an  initial  angle  of  attack,  the  development  of  a 
shock  wave  at  the  leading  edge  will  affect  the  aerodynamic  coefficients.  The 
magnitude  of  this  effect  is  demonstrated  by  evaluating  the  aerodynamic 
coefficients  Lj^,  M^^,  L^,  for  an  oscillating  flat-plate  airfoil  at  a  positive 
initial  angle  of  attack.  The  calculations  are  performed  for  a  free -stream 
Mach  number  of  5  and  for  M  0  equal  to  0.5,  0.9,  1.3,  and  1.7.  The  results 

GO 

of  the  calculation  are  shown  in  Figs.  2a,  2b,  and  2c.  The  importance  of  this 
nonlinear  aerodynamic  effect  is  immediately  apparent  from  these  figures. 

The  effect  is  significant  even  at  values  of  M^O  as  low  as  0.5. 


Numerical  Examide 


0 


1.0 


t.o 


Fig.  2b.  Oscillatory  aerodynamic  coefficients  (continued). 
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0  1.0  2.0 

M.2 


Fig.  2c.  Oscillatory  aerodynamic  coefficients  (continued). 


V.  OSCILLATORY  AERODYNAMIC  INFLUENCE  COEFFICIENTS 


In  certain  instances  it  may  be  convenient  to  describe  the  motion  of  the 
airfoil  strip  in  terms  of  displacements  rather  than  in  terms  of  angles  and 
displacements.  The  aerodynamic  lift  L,  leading  edge  moment  M.  and  hiimte 
moment  T  are  replaced  by  an  equivalent  system  of  aerodynamic  forces.  The 
displacements  {h}  and  equivalent  aerodynamic  forces  {f}  are  related  by 
means  of  oscillatory  aerodynamic  influence  coefficients.  The  matrix  [  Cj^l 
of  aerodynamic  influence  coefficients  (AICs)  is  defined  by  the  equation 

{f|  =  p^«V8[Cj^]lh|  .  (29) 

where  b^  is  the  reference  semichord  and  s  is  the  wing  semispan.  The  aero¬ 
dynamic  influence  coefficients  can  be  obtained  directly  from  the  oscillatory 
aerodynamic  coefficients. 

Consider  the  equivalent  force  system  shown  in  Fig.  3b.  The  equivalent 
forces  are  arbitrarily  placed  at  the  quarter-chord,  control  surface  hinge  line, 
and  trailing  edge.  The  airfoil  strip  is  assumed  to  be  of  unit  width,  and  the 
aerodynamic  coefficients  are  defined  as  in  the  preceding  section.  The 
relationship  between  the  matrix  of  AICs  and  the  aerodynamic  coefficients  can 

4 

be  shown  to  be 


(1  +  b/2d) 

-b/d 

(b/c^)(3b/2d  -  1) 

[Chl 

=  (4/s)(b/b^)2 

-b/2d 

b/d 

-(b/c^)(3b/2d) 

0 

0 

b/Ca 

V 


(1  +  b/2d) 
-b/d 
b/d 


-b/2d 

b/d 


-(b/d  +  b/Cj^) 


0 

0 

b/c. 
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where  denotes  the  control  surface  chord  and  d  is  the  distance  between 
the  forward  and  aft  control  points. 


Fig.  3.  Original  (a)  and  equivalent  (b)  force  systems 
and  geometry  for  the  oscillatory  case. 


APPENDIX  A 


A  Particular  Solution  of  the  Hypersonic  Small  Di»turbance 
Equatiom  for  Two-Dimen>ional  Un»t«*dy  Flow 

Consider  the  small  transverse  motions  of  a  two-dimensional,  slender, 
pointed  body  immersed  in  a  steady,  uniform  stream  of  a  perfect  gas 
(see  Fig.  1).  The  fluid  is  inviscid  and  does  not  conduct  heat.  An  order  of 
magnitude  analysis  of  the  equations  governing  the  disturbances  introduced  by 
the  presence  of  the  body  may  be  carried  out  under  the  assumption  that 

M  »  1  ,  T  «  I  ,  Mt~  0(1)  ,  k*  ~  0(1)  .  (A-1) 

where  t  is  a  measure  of  the  slope  of  the  body  surface  and  where  M  denotes 
the  Mach  number.  The  parameter  k*  =  L/U^T*,  where  T*  is  a  character¬ 
istic  time  associated  with  the  time -dependent  transverse  motions  and  L  is  a 
characteristic  length  in  the  streamwise  direction,  is  a  measure  of  the 
unsteadiness  of  the  flow.  Introduction  of  the  appropriate  variables^  into  the 
equations  of  motion  and  neglect  of  the  higher  order  terms  leads  to  the  field 
equations  of  the  hypersonic  small  disturbance  theory; 


p.  +  U  p  +  (pv)  =  0 
t  oo'^x  '  'y 

(A- 2a) 

p(v.  +  U^v  +  vv  )  +  p  =0 
t  oo  X  y'  ’^y 

(A -2b) 

S*+U  S  +S  =  0 
t  oo  X  y 

(A.2c) 

p(u.  +  U  u  +uu  +  vu)+p  =  0  , 

'  t  GO  X  X  y'  '^x 

(A -2d) 
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Mtdiere  S  =  p/p^.  The  preeaure  end  density  are  denoted  by  p  and  p, 
respectively.  The  perturbation  velocities  along  the  x  and  y  axes  are  u  and 
V.  The  free  stream  velocity  is  denoted  by  U^. 

The  boundary  condition  on  the  body  is 

v(x.  0.  t)  =  +  Yj  =  U^6(x.  t)  .  {A-3) 

where  Y(x,  t)  defines  the  body  surface.  The  disturbances  are  assumed  to 
vanish  far  ahead  of  the  body.  The  changes  in  the  flow  conditions  across 
shock  waves  are  governed  by  the  appropriate  jump  conditions  (see  Ref.  5). 

Equations  (A-2a),  (A-2b).  and  (A>2c)  can  be  solved  independent  of 
Eq.  (A -2d),  and  these  equations,  together  with  the  appropriate  boundary 
conditions,  constitute  the  reduced  hypersonic  problem.  A  steady  flow 
problem,  equivalent  to  this  reduced  hypersonic  problem,  may  be  constructed 
by  means  of  the  following  transformation: 


•  ♦  * 

■  *  'TJ—  • 

(A-4a) 

oo 

X  =  X  , 

(A-4b) 

7=  y  . 

(A-4c) 

The  equations  governing  the  equivalent  steady  problem  are 


(A-5a) 

(A-Sb) 


.28. 


(A-5c) 


where  3f  =  p/?"^; 

|f(X;T)  =  p(x,  t) 
p'(3E:T)  =  p(x,  t) 
v(x:T)  =  v(x,  t) 

The  boundary  condition  on  the  body  becomes 

^(3E;T)=  U^?’(Sr:r) 


Assuming  that  the  flow  processes  are  isentropic  in  shock-free  regions, 
may  combine  Eqs.  (A-5)  into  the  following  form: 


[^Qo  ^  ^  ° 


where 

+  2S'/(y  -  1), 

e  =  zs/iy  - 1), 

a  =  -YP/P- 


(A -6a) 
(A -6b) 
(A-6c) 

(A -7) 

we 

(A-8a) 

(A-8b) 
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Considering  the  upper  surface  of  the  body,  we  employ  the  following  solution  of 
Eqs.  (A-8)  for  a  region  r  of  isentropic  flow. 

F  =  constant  along  curves  defined  by  d^:  dSf  =  (v  -f  a^: 

G  =  constant  throughout  the  region  T. 

This  choice  of  solution  implies  that  the  disturbances  reflected  back  to  the 
body,  which  arise  due  to  the  presence  of  curved  shock  waves,  are  negligible. 

The  pressure  distribution  in  F  is,  therefore, 

~  2y/(y-1) 

(  V  -  V  \  ' 

l+^S - a^)  ' 

n  ' 

where  v  ,  ff  ,  S’  are  the  conditions  at  some  reference  point  in  the  region  F. 
n  *^n  n  r  a 

The  pressure  distribution  on  that  part  of  the  body  surface  contained  within 
F  is 

1  +X^M^(6  .  6^)]  (A- 10) 

If  F  is  the  region  downstream  of  the  bow  shock  (assumed  to  be  attached  to  the 
nose)  and  upstream  of  any  successive  shock  waves  emanating  from  the  body, 
then  the  reference  point  may  be  taken  at  the  nose  of  the  body  and  the  conditions 
M^,  p^  determined  from  oblique  shock  theory.  The  solution  for  such  a  case, 
when  transformed  back  to  the  x,  y,  t  coordinate  system,  is 

2y/(y-1) 

p(x.  t)  =  p^(t*)  1 1  +  M^(t*)l6(x,  t)  -  6n(t*)])  (A- 1 1) 


5 

\ 


The  quantities  denote  the  aerodynamic  pressure, 

Mach  number,  and  flow  deflection  angle  at  the  nose,  respectively,  evaluated 

at  the  time  t*  =  t  -  x/U„.  In  the  case  that  there  is  no  shock  wave  at  the 

CO 

nose,  Eq.  (A- 11)  simplifies,  and 

I  v-1  iZy/Cy-i) 

P(x.  t)  =  Pqo[1  + 

The  solution  (A- 11)  may  be  considered  as  the  unsteady  version  of  the  shock 
expansion  formula  for  the  case  of  two-dimensional  hypersonic  flows  over 
slender,  pointed  bodies  a  small  angles  of  attack.  To  recapitulate  briefly, 
this  result  was  developed  under  the  assumptions  that 

1.  The  fluid  is  a  perfect  gas, 

2.  The  effects  of  fluid  viscosity  and  thermal  conductivity  are 
negligible, 

3.  The  flow  is  such  that 

M  »  1  ,  T  «  i  ,  Mt  ^  0(1)  ,  k*  --  0(1) 

4.  The  fluid  processes  in  shock-free  regions  are  isentropic, 

5.  Disturbances  reflected  back  to  the  body  are  negligible. 
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APPENDIX  B 


Definition  of  Coefficients 


Pressure  Distribution  on  an  Airfoil  Stric 


The  coefficients  appearing  in  Eqs.  (19)  tiirough  (22)  are  defined 


as  follows: 


Poo^l 

Poo^2  =*'2u+*l2i 
Poo^3  ■^'lu'^^Ui 
PC0G4  =*l2u  +  »\l 
Pqo^S  *'51 
Poo®6  “  *'61 
PcoS 

PqoS  "  *l7u  *'2i 
Poo^9 

Poo^l0**''6u 
Pencil  =  3ll3u  +  *  If 
Pqo^IZ  "  *'21 
PaoGl3  =  *'3u  +  *i4f 
Poo^l4  ^  *71 
Pa.Gi5  =  *W  +  *!ll 


Poo°l6  ''*9u 


*  lOu 

The  additional  subacripts  u  and  i  indicate  that  quantities  are  to  be  evaluated 
on  the  upper  and  lower  surfaces,  respectively. 

The  coefficients  appearing  in  Eqs.  (26)  and  (27)  are  defined  as 
follows: 

Hj^(m,n)  =  ik  +  (ik  + 


H2(m,„)=  (l  +  ik|)0„t[l  -Ik  ]g„  , 

Hjte,,.  r.B)  =  ik  Gp  t(lk  |)g^  k(ik  l^j(G^  +  G.) 


2  2 

r,  .)  .{1  t  Ik  «)Gp  k  [l  -  ik  I  -  ^(1)  ]g^ 

k[(l  t  ik  k(i^  -  ik)|  -  +  °.l  • 

H^ip.  r)  =  t  ik  k  G^  jl  t  ik  k  5(k2  ^  -  u) 


k 

T 


f] 
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Oscillatory  Aerodynamic  Coefficients 


The  coefficients  appearing  in  Eqs.  (28)  are  defined  as  follows: 


bK^(m.  n)s  J  Hj(m,  n)dx  , 

J^2b 

'  Hj(p,  q,  r,  s)dx  , 
bK,(m,  n)  =  /  H,(m,  n)dx  , 

Jo 

J-2b 

f  H.(p,  q.  r.  s)dx 
’St 

/.2b 

bKg(p,  r)  =  /  HgCp,  r)dx 
•'*b 

2 

b  K/(m,n)  =  /  xH,(m,n)dx  , 

^  Jo  ‘ 

2  ^2b 

b  Ky(p,q,  r,8)=  I  xH2(p,  q,  r,  8)dx 

J^ 

2 

b  KQ(m,n)  =  #  xH,(m,n)dx  , 

^  Jo  ^ 

2  ^  2b 

b  KQ(p,q,  r,s)=  I  xH.(p,  q,  r,  8)dx 

J\ 

b  KjQ(p,  r)=y  xHgtp,  r)dx 
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